SUSY contributions to the charge 
asymmetry in K ± — ► , K ± t r t~ decays 
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! We analyse the contributions to the charge asymmetry in — > ir £ + £~~ decays induced 
by gluino-exchange diagrams in the context of supersymmetric models with generic flavour 
couplings. We show that sizeable deviations with respect to the Standard Model are 
possible only under special circumstances. Within this scenario we set an upper limit of 
O-i. about 10~ 3 for the relative charge asymmetry - integrated for M £ +£- > 2M n - of both 
muon and electron channels. We also show that this limit is close to saturating a model- 
independent upper bound on the charge asymmetry derived from the present constraints 
on T(K L - jr»< 



I. CP violation is one of most interesting and least known aspects of particle physics. The 
recent results of KTeV [|J and NA48 [fj] about e' / e have unambiguously shown the existence 
of CP violation in \AS\ = 1 transitions (the so-called direct CP violation) and ruled out 
superweak scenarios. The experimental measurements of e'/e are generally compatible with 
the theoretical expectations within the Standard Model (SM) ||. Nevertheless, the latter 
are affected by sizeable theoretical uncertainties, and it is difficult to constrain non-standard 
effects. At the moment the theoretical uncertainties on e'/e are so large that we cannot even 
exclude that this observable is completely dominated by new physics (NP) contributions. 
Given this situation, it is highly desirable to obtain new independent information about 
CP violation in |AS| = 1 transitions. 

Among direct-CP-violating observables, particularly interesting are those accessible in 
rare decays. In these processes the smallness of the SM contribution leads to identifying 
more easily a possible large NP effect, whereas the simplicity of the hadronic structure helps 
us to keep under control the theoretical uncertainties [|J . An observable that satisfies these 
requirements is the charge asymmetry in — > 7r ± £ + £, defined as 



where q 2 = M£ +i - is the dilepton invariant mass and I = e or fi. 

This asymmetry is a pure direct-CP-violating observable. A non-zero A^ is generated 
by the interference between the absorptive contribution of the long-distance amplitude and 
a CP- violating phase of short- distance origin ||. As we shall discuss below, the kinematical 
cut on the dilepton invariant mass (q 2 > Sq) is a useful tool to maximize the CP-violating 
effect; indeed, the SM expectations of this asymmetry are given by || 
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The smallness of these figures leads to considering this observable as a good probe of 
possible NP effects. One of the most promising scenarios where sizeable non-standard CP- 
violating contributions could be generated is the supersymmetric (SUSY) extension of the 
SM with generic flavour couplings and minimal particle content. In particular, it has been 
recognized that in this context gluino-mediated penguin diagrams could naturally account 
for the observed value of e'/e Q . The purpose of this paper is to analyse the consequences 
that this scenario could have on A^. We shall therefore assume that the CP-violating 
phase of the s — ► d£ + £~ amplitude is entirely dominated by gluino-mediated penguin 
diagrams or, to be more specific, by the contributions of the dimension-5 chromomagnetic 
and electromagnetic dipole operators (CMO and EMO). The couplings of these operators, 
determined by the mismatch between quarks and squarks mass matrices, appear also in e 
and e' . We will therefore extract the allowed range of these couplings from the measured 
values of e and e' in order to analyse the possible effects on A^. 

Our conclusions are that, in the general case, SUSY effects are at most as large as 
the SM results in (§) and thus not particularly interesting. Only assuming a cancellation 
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between two independent SUSY contributions to e', or in a fine-tuned scenario, is it possible 
to reach higher values. In any case, within the minimal SUSY extension of SM considered 
here, cannot exceed the 10~ 3 level (independently of the q 2 cut). We will also show that 
this bound is very close to a model-independent limit on the charge asymmetry which was 
extracted, by means of isospin symmetry, from the present constraints on T(Kl — > tt°£ + £~). 

This paper is organized as follows: we shall first discuss the amplitude decomposition 
of K — > n £ + £~ decays and the generic expression of the charge asymmetry; then we 
will introduce the effective Hamiltonian describing the SUSY short-distance contributions; 
finally we shall discuss the phenomenological constraints on the SUSY phases and the 
corresponding bounds for the charge asymmetry. 



II. The charge asymmetry is produced by interference between the CP-conserving strong 
phase and the CP-violating weak phase of the decay amplitude. The latter is generated 
by the exchange of heavy particles (short-distance), the former is due to the K — > 3tt 
intermediate state and thus belongs to the long-distance part of the amplitude (K ± — > 
37r — ► 7r ± £ + £ _ ). As shown in ||, the decay width of K ± — > tt £ + £~ is largely dominated 
by long-distance effects and, as long as CP violation is not considered, short-distance 
contributions can be neglected. 

In the long-distance part of the K — > ir £ + £~ amplitude the two leptons are always 
produced by a virtual photon (K ± — > 7T7* — > 7r ± £ + £ _ ). Since we are interested in the 
interference between short- and long-distance terms, all short- distance contributions where 
the two leptons are not in a vector state can be neglected. We can therefore parametrize 
the decay amplitude in terms of a single vector form factor, W(z), defined as || 

A(K ± (k)^7r ± (p)e + (i + )£(n) = -^w(z)(k+ P ru(i + ) lfi v(n , w 

where q = k — p and z = q 2 jM\. Integrating this amplitude over the phase space leads to 

where 4r| < z < (1 — r^) 2 , = mi/M K and p(z) = Jl — Ar\l z{\ + 2r 2 jz)}?l 2 {\, z, r 2 ). 
Referring to the discussion in Ref . 0] , we decompose the form factor as 

W {z) = W*>° l {z) + y^n) • (5) 
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Here PU 7r7r (z) denotes a non-analytic contribution, largely dominated by the dipion inter- 
mediate state, as a consequence of the small q 2 of the lepton pair, which exhibits a branch 
cut along the real axis starting at z — 4r 2 _. The polynomial term W pol (z) includes both 
long- and short-distance contributions. The latter are completely dominant in the real 
part of W pol (z), which can be fitted from experiments. Parametrizing the polynomial term 
as W pol (z) = GpM\[a + bz), consistently with a chiral expansion at the next-to- leading 
order ||, the values of a and b fitted by the BNL E865 Collaboration are reported in Ta- 
ble [I]. On the other hand, the imaginary part of W po1 can be produced only by the direct 
CP-violating weak phase of short- distance origin. 
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(-0.587 ±0.010) 


(-0.655 ± 0.044) 


(-20.6 ±0.5) x 10~ 8 


(-2.4 ± 1.2) x 10~ 8 



Table 1: Experimental values of the parameters determining W po1 ||/ and W 7rw (z) j§, 0/ 



The function W n7T (z) has already been computed in the literature |5], [| and we report 
it here for completeness: 
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where the values of a and (3 (fitted from K — > 3ti decays) are shown in Table[l], z Q = 1/3+r 2 , 
$(z) — 1 ± z/ry (r v ~ 2.5) and x{z) is defined byQ 
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The interference between the CP-violating phase in W /po ' and the absorptive contribu- 
tion to W(z) n7T leads to the following difference between the widths of the charge-conjugated 
modes 



rt-r: 
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Obviously, in order to obtain an adimensional asymmetry we should normalize the width 
difference to the sum of the widths. However, we stress here the strong dependence of 
this observable from the dilepton invariant mass: owing to the kinematical threshold in 
the absorptive part, only when the two-pion intermediate state can be on shell is the 
asymmetry different from zero. This observation suggests the construction of integrated 
asymmetry setting a cut on q 2 above 4mJ, i.e. setting so = 4mjj. in Eq. (0). Since electron 
and muon channels have a very similar phase space for q 2 > Am 2 , with such a cut the 
normalized asymmetry turns out to be very similar in the two modes. On the other hand, 
since most of the available phase space for the electron channel is below this cut, from the 
experimental point of view the muon channel appears a better candidate for the study of 
this asymmetry. 



III. As discussed in the previous section, we look for short- distance contributions to the 
s — > d£ + £~ amplitude, with the lepton pair in a vector state and possibly with a sizeable 
new weak phase. As shown in Ref. flOfl , in the presence of non-minimal flavour mixing, 



1 x satisfies the relation: x(0) = 0; this condition follows from the possibility to include x(0) into a 
redefinition of the polynomial form factor. 
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the largest SUSY contributions to this type of amplitude are produced by the dimension- 
5 CMO and EMO. All the other contributions are in fact naturally suppressed once the 
bounds from other processes are taken into account (barring accidental cancellations). The 
structure of the effective Hamiltonian necessary to describe these contributions reads 
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with the operators expressed in the following basis 
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Note that in this basis the operators have well defined properties under parity transforma- 
tions: Qf induce parity-conserving transitions and parity-violating ones. 

The Wilson coefficients of these operators induced by gluino exchange are given by 
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Here (Slr)^ = (M 2 ) iL j R /m c j denotes the off-diagonal entries of the down-type matrix in 
the super-CKM basis, and x = mhjmh the ratio of gluino and (average) squark masses 



squared. The explicit expression of the loop functions F(x), G(x) can be found in |10 



We have also considered the CMO, even if it does not participate at tree level, because of 
the large mixing between CMO and EMO induced by QCD interactions at the one-loop 
level. Taking into account the 2x2 anomalous-dimension matrix of these operators |T0f , 
computed at lowest order, the Wilson coefficients evolved down to charm scales reads 



?? 2 [C 7 K) + 8(1 -rT^K)], 
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Starting from Eq. 

as 



where 



<\An>i)J \a s {m b ) ) \a s (m c ) J " \a s (500 GeV) 

151), we find it convenient to rewrite the Wilson coefficient of the EMO 
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(5 LR ) sd ± (5 LR )* ds . 
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In order to identify the SUSY contribution to the polynomial form factor, we need to 
evaluate the EMO matrix element between K + and tt + external states: 



{u + \Q+\K + ) = 2% 



16vr 2 M K 



Pfl (n + )p u (K + )F^. 



(21) 



As is usually done in the literature, we have expressed the hadronic matrix element in 
terms of a suitable parameter [|H], [12]], expected to be 0(1), which encodes the non- 



perturbative dynamics. This parameter has recently been computed on the lattice ]14 



confirming the estimate ~ 1 made in Ref . |T2[ . We shall leave it as a parameter through 
all the paper, but for completeness we report here the recent lattice result 



B T (n = 2 GeV) = 1.21 ± 0.09 ± 0.04^ 



0.07 
00 



[14]. 
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Now using Eqs. (|18|)-(|21~|) we can write the full matrix element of the Hamiltonian ( |I~0| ) 
relevant to K ± —>■ n ± £ + £~ decays 
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Then, according to the definition of W(z) in Eq. 
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\Y (1)J V m 
g — 500 GeV) = 0.39. Inserting this result in Eq. 
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(25) 
we can 



1.0 ±0.1) 



'Y(x)\ /500 GeV\ / a s (m § ) 
Y (l) ) \ fit ) V« s (500 GeV) 



xB T \Im(S + )\, (26) 



where the error includes the uncertainty in the experimental parameters. 



IV. Before starting a numerical analysis of the SUSY contribution to the charge asym- 
metry, we discuss here a more general upper bound on |Im(IU poZ )| (and thus on the charge 
asymmetry), which can be extracted from the experimental upper bound on T(Kl — > 
^ £+£-). 

Isospin symmetry relates in a model-independent way the short-distance components 
of — ► tt ± £ + £~ and —>■ ix ± £ Jr £~ amplitudes. Indeed, in both cases, the hadronic 
current is necessarily a AI = 1/2 operator of the type sTd (or dTs). In the Kl case, 
the approximate CP-odd combination of K° and K° states select the imaginary part of 
this amplitude, which can therefore be constrained by using the experimental bounds 
on T(Kl — > 7T £ + £~). Assuming that long-distance contributions to Ki ix ± t r £~ are 



5 





SM 


SUSY 


Mo del- Independent 


|A e (4m;f)| 


0.07 x |Im(A t )| 


0.13 x |Im(5+)| 


1.3 x 10" 4 


|A e (4m*)| 


0.53 x |Im(A t )| 


1.0 x \lm(5 + )\ 


1.0 x 10" 3 


|A,(4m 2 J| 


0.21 x |Im(At)| 


0.40 x |Im(5+)| 


6.4 x 10" 4 


|A„(4m*)| 


0.54 x |Im(Ai)| 


1.0 x |Im(5+)l 


1.6 x 10-3 



Table 2: Summary of the results for both the electron and muon channel asymmetry, 
evaluated in the SM ^ and in the MSSM, with and without the kinematical cut on the 
dilepton square mass. 



negligible (which is certainly a good approximation if T(Kl — > ~ 
current experimental bound Q), we can write 
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where tl is the mean lifetime of the Ki and the integral on p(z) extends to the whole 
phase space. Using this relation we obtain the following model-independent upper limits 
on the charge asymmetries of electron and muon modesQ: 



A e (4m^ 



< 43 x y/BR(K L 

< 82 x ^BR(K L 



< 1.0 x 10' 



< 1.6 x 10" 



(28) 
(29) 



where the numerical values have been obtained using the experimental bounds BR(Kl — > 
7r°e + e') < 5.8 x lO" 10 and BR(K L -> ir°u, + u,-) < 3.8 x 10~ 10 obtained by the KTeV 
Collaboration PB| . 

In Table [2| we compare the model-independent upper bounds for the charge asymmetries 
with SM and MSSM expectations, expressed in terms of the imaginary parts of the CKM 
factor X t = V t * s V tc i and the SUSY parameter 6 + , respectively. Since Im(A t ) ~ 10~ 4 , it is 
clear that the SM could only account for a few per cent of the general bounds. From this 
perspective it is particularly interesting to understand how large the SUSY contribution 
could be. 

The bounds on lm(5 + ) [and more in general on the lm(5LR)i/s\ have been widely 



discussed in the recent literature (for extensive discussions, see Ref. Jl5| and references 
therein). In general there are two classes of bounds: those coming from |e| and BR(Kl — > 
7r°e + e~) (direct bounds) and those extracted from Re(e'/e) (indirect bounds). In the 
first case lm(S + ) is directly involved, whereas in the second case - dealing with parity- 
violating transitions - only lm(S~) is directly involved. As clearly shown by Eq. (p0|), 
6 + and S~ are naturally related: if we assume these two quantities to be of the same 
order of magnitude, then the measurement of Re(e'/e) leads to a bound on |Im(<5 + )| of 

2 We report separately the bounds on electron and muon modes, obtained by the corresponding Kl — * 
tt°£ + £~ decay, in order to take into account possible violations of lepton universality. 
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Figure 1: Electron differential asymmetry as a function of z. 



about 10 -6 0. In this case the SUSY contribution to the charge asymmetries could be 
at most of the same order as the SM one. On the other hand, if we allow a fine-tuned 
scenario where |Im(5 + )| 3> |Im(5~)|, then only the (weaker) direct limits from |e| and 
BR(Kl — > 7r°e + e~) have to be satisfied. The constraints from BR(Kl — > 7r°e + e~) are 
equivalent to the model-independent bounds reported in Table |2] and thus they certainly 
allow for sizeable deviations with respect to the SM. In the case of |e| there is a considerable 
uncertainty, induced by hadronic matrix elements: employing optimistic values for the 
latter (or allowing cancellations with other SUSY effects) it is possible to obtain bounds 
on Im(5 + ) above 10~ 4 [15|. Under these special circumstances the SUSY contribution to 



the charge asymmetry could become substantially larger than the SM one, reaching values 
close to the model-independent bounds. 

Before concluding, we point out that an observable particularly useful from the exper- 
imental point is the differential asymmetry, defined by: 



6AAz) 



dT^(z)/dz — dT e (z)/dz 
dT^ (z) jdz + dTj (z) / dz 



(30) 



This observable could be particularly useful in an experimental set up where it is difficult 
to obtain a precise flux normalization and thus precise width measurements. In Figure [l] 
we plot SAi(z), normalized to the SUSY phase Im(<5 + ), as a function of z. 



V. In this letter we have analysed the charge asymmetry in K ± — > ix ± t r t~ decays in the 
framework of supersymmetric models with minimal particle content and generic flavour 
couplings. We have shown that in general the expectations are of the same order of 
magnitude as in the SM case. Nevertheless, under special circumstances it is possible to 
relax the indirect constraints on the SUSY CP-violating phases and obtain results in the 
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range 10~ 4 < A^(4m^.) < 10~ 3 , or above the SM expectation. We have also shown that 
under general assumptions A^Arn^) cannot exceed the 10~ 3 level. This model-independent 
bound is derived, using isospin invariance, from the experimental constraints on T{Kl — > 



7T° 



The first measurement of the charge asymmetry in the muon mode has recently been 
announced by the HyperCP collaboration [16]: 

A M (4mJ) = -0.02 ± 0.11 (stat) ± 0.04 (syst) . (31) 

Unfortunately at the moment the sensitivity is very far from the interesting region, still well 
above the model-independent bound. In view of possible improvements of this measure, 
we can summarize as follows the possible outcome of a non- vanishing result: 

• A^(4m^) ^ 10~ 4 : we could conclude that the charged kaon system is well described 
by the SM. Moreover, we could put a new strong constraint on the SUSY phase 
lm(5 + ). 

• 10~ 4 < A^(4m^.) < 10~ 3 : we would have a clear signal of new physics, compatible 
with the results obtained within the Minimal Supersymmetric Model with generic 
flavour couplings. 
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